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Abstract 

Bosons with non-zero spin exhibit a rich variety of superfluid and insulating phases. Most phases 
support coherent spin oscillations, which have been the focus of numerous recent experiments. 
These spin oscillations are Rabi oscillations between discrete levels deep in the insulator, while deep 
in the superfluid they can be oscillations in the orientation of a spinful condensate. We describe 
the evolution of spin oscillations across the superfiuid-insulator quantum phase transition. For 
transitions with an order parameter carrying spin, the damping of such oscillations is determined 
by the scaling dimension of the composite spin operator. For transitions with a spinless order 
parameter and gapped spin excitations, we demonstrate that the damping is determined by an 
associated quantum, impurity problem of a localized spin excitation interacting with the bulk critical 
modes. We present a renormalization group analysis of the quantum impurity problem, and discuss 
the relationship of our results to experiments on ultracold atoms in optical lattices. 

PACS numbers: 05.30.-d, 03.75.Kk, 71.10.-w 
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I. INTRODUCTION 



An important frontier opened by the study of ultracold atoms has been the investigation 
of Bose-Einstein condensates of atoms carrying a nonzero total spin F. The condensate 
wavefunction then has rich possibihties for interesting structure in spin space, analogous 
to structure investigated earlier in superfluid ^He. The dynamics of the atomic conden- 
sate is described by a multi-component Gross-Pitaevski (GP) equation, and this allows for 
interesting oscillations in the orientation of the condensate in spin space. Such coherent 
spin oscillations have been observed in a number of recent studies of F = 1 and F = 2 
condensatesii^'^i^i^i^. 

Coherent spin oscillations of a rather different nature are observed in the presence of a 
strong optical lattice potential^i^. Here there is no Bose-Einstein condensate and the ground 
state is a Mott insulator: each minimum of the optical lattice potential traps a fixed number 
of atoms, and tunneling between neighboring potential minima can be ignored. Now the spin 
quantum number leads to a number of discrete atomic levels associated with each minimum, 
with the energies determined by the "cold collisional" interactions between the atoms. The 
coherent spin oscillations are then the Rabi oscillations between these atomic levels^i^. 

Note that the spin oscillations in the Mott insulator emerge from a solution of the full 
quantum Schrodinger equation in the finite Hilbert space of each potential minimum. In 
contrast, the oscillations in the superfluid condensatesi'^i'^i^'^ are described by classical 
GP equations of motion obeyed by the multicomponent order parameter, representing the 
collective evolution of a macroscopic condensate of atoms. 

It is the purpose of this paper to connect these distinct spin oscillations to each other 
across the superfiuid-insulator quantum phase transition. The equilibrium properties of the 
superfiuid-insulator transition of spinful bosons are quite complicated, and the very rich 
phase diagram has been explored theoretically, both for— ^iiii^ F = 1 and for— F = 2. It is 
not our purpose here to shed further light on the nature of this phase diagram, or on the 
possibilities for the experimental realization of the various phases. Rather, we will examine 
representative cases which display the distinct possibilities in the evolution of coherent spin 
oscillations. 

Deep in the superfluid, the spin oscillations are associated with normal modes of the GP 
equations, about the superfluid state. Some of these modes carry spin, and so will contribute 
to the spin oscillations. If the superfluid state breaks spin rotation invariance, then such 
a spin-carrying mode will be gapless. Otherwise, the spin oscillations are gapped, ie, they 
occur at a finite frequency. The dominant damping of these oscillations will arise from the 
creation of low energy 'phonon' excitations of the superfluid. However, the gradient coupling 
to such 'Goldstone' modes will suppress the decay, and so one expects the oscillations to be 
well-defined. 

Deep in the Mott insulator, the discrete energy levels of each potential minimum lead to 
undamped oscillations. Coupling between these minima is expected to again lead to weak 
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damping. 

Assuming a second-order quantum phase transition between these hmits, we can expect 
that the spin oscillations will experience a maximum in their damping at the critical point, 
and possibly even cease to exist as well-defined modes. There are a plethora of low energy 
excitations at the quantum critical point, and their coupling to the spin modes is not con- 
strained by the Goldstone theorem. The primary purpose of this paper is to describe this 
enhanced damping in the vicinity of the quantum critical point. 

Our analysis shows that behavior of the spin oscillations falls into two broad classes, and 
we will present a detailed analysis of a representative example from both classes. The classes 
are: 

(A) The order parameter for the quantum transition carries non-zero spin. The spin excita- 
tion spectrum is gapless at the quantum critical point, and the spin operator is characterized 
by its scaling dimension. Typically, the spin operator is a composite of the order parameter, 
and standard methods can be used to determine its scaling dimension. The value of this 
scaling dimension will determine the long-time decay of spin correlations, and hence is a 
measure of the damping. 

(B) The order parameter of the superfluid-insulator transition is spinless. In this case, it 
is likely that all excitations with non-zero spin remain gapped across the quantum critical 
point. We then have to consider the interaction of a single gapped spin excitation with the 
gapless, spinless critical modes. Such a problem was first considered in Refs. [l^EsI for the 
case of gapped fermionic excitations. Here we will show that closely related considerations 
apply also to the present case with gapped bosonic excitations. The dispersion of the gapped 
excitation is argued to be an irrelevant perturbation, and so to leading order one need only 
consider the coupling of a localized spin excitation interacting with the bulk critical modes. 
This gives the problem the character of a quantum impurity problem. Depending upon 
whether the localized-bulk coupling is relevant or irrelevant, we then have two sub-cases. 
(Bl) If the coupling between the localized and bulk excitations is relevant, then a renormal- 
ization group (RG) analysis is necessary to understand the structure of spin correlations. A 
brief account of this RG analysis was presented earlier—, and here we will present further 
details describing the new, non-trivial impurity fixed point controls the long-time physics. 
(B2) For the case of an irrelevant localized-bulk coupling (which we will also present here), 
the damping is controlled by scaling dimensions of the bulk theory, and no new impurity 
dimensions are needed. 

The representative examples we will consider in this paper will all be drawn from the 
simplest case of F = 1 bosons in an optical lattice potential, with an even number, A^, of 
bosons per site^^iiiii^. The mean-field phase diagram of a model Hamiltonian for F = 1 
bosons is shown in Fig. [1], along with the locations of the quantum phase transitions of the 
two classes described above. Further details on the phases and phase transitions appear in 
the body of the paper. We expect that transitions with the more complex order parameter 
categories possible for— F = 2 and higher—, will also fall into one of the categories we have 
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FIG. 1: Phase diagram for the Hamiltonian in Eq. (j2.ip calculated using the mean-field theory of 
Section |TT1 The three phases that are included are the spin-singlet insulator (SSI), the spin-singlet 
condensate (SSC), and the polar condensate (PC). The arrows labeled 'A' and 'B' denote the two 
classes of transition defined in SectionlH The calculation has been performed using quantum rotors, 
corresponding to the canonical ensemble with the filling factor an even integer. The horizontal and 
vertical axes give the tunneling strength t and the spin-dependent interaction J, both in units of 
the spin-independent part of the interaction, U. 

described above. 

We will begin in Section [TTl by describing the model Hamiltonian and a simple mean- field 
theory that can be used to treat those phases with which we shall be concerned. We defer a 
more thorough description of the symmetries and possible phases of the model to Section HTH 
where we consider the continuum limit of the theory. 

In Section IIVI we use this continuum theory to describe the low-energy properties of the 
phases with which we are concerned, focusing on the behavior of spin excitations. We then 
turn, in Section |Vl to the behavior across transitions falling into the classes identified above. 
The remaining sections treat these cases in more detail. 



II. MODEL AND MEAN-FIELD THEORY 
A. Model 

Bosonic atoms trapped in an optical lattice potential, at sufficiently low temperatures that 
all atoms occupy the lowest Bloch band, are well described by the Bose-Hubbard model. 



(See, for instance, Ref. Il8| for a recent review.) The extension of this model to the case 
where the bosons carry spin is straightforward. While most of our results apply for general 
spin F 7^ 0, we will treat explicitly the case F = 1 and draw attention to the generalizations 
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when appropriate. The Hamiltonian can then be written as 



(2.1) 



where summation over repeated spin indices /i is imphed throughout. The operator 
annihilates a boson at site i with spin index fi G {x,y,z}] this basis will be the most 
convenient for our purposes. 

The first term in Ti, involves a sum over nearest-neighbor pairs of sites, and allows the 
bosons to tunnel, with 'hopping' strength t. We will restrict attention to square and cubic 
lattices, but the results can be straightforwardly generalized to other cases. The function 
V{n) contains an on-site spin-independent interaction and the chemical potential, and can 
be written in the form 



In the final term of Eq. (12. ip . Lj is the total angular momentum on the site i, given by 
(for F = 1) Li^p = —ie^upblpbiu, where e^^p is the completely antisymmetric tensor. For 
F = 1, this term is the most general quartic on-site spin-dependent interaction. For spin 
F, the boson operator becomes a tensor of rank F and there are F + 1 independent quartic 
interaction terms corresponding to different contractions of the spin indices. We will not 
include any direct interactions between spins on neighboring sites, nor long-range polar 
forces between the atoms. 

Suppose that V{n) has its minimum near some even integer, N, and the couplings are 
tuned so that the model is particle-hole symmetric around this filling. Requiring this sym- 
metry corresponds to restricting consideration to the case of integer filling factor, and is 
equivalent to using the canonical ensemble^^. 

The spinless Bose-Hubbard model, with a single species of boson, has a transitioni^i^ 
from a Mott insulator when U ^ t to a, superfluid when U <^ t. When the bosons have 
spin, various types of spin ordering are possible within both the insulator and superfluid. 

With an even number of particles per site, the simplest insulating phase, the spin-singlet 
insulator (SSI), has a spin singlet on each site. This will be favored energetically when 
J > 0, and we will concentrate on this case in the following. For J < or odd A^, the net 
spin on each site will be nonzero, and the system will be well described by a quantum spin 
model, allowing various forms of spin ordering within the lattice^. 

In a simple superfluid, the bosons condense, so that {bi^) ^ 0, breaking both gauge and 
spin symmetries. For J > 0, a so-called polar condensate (PC) is favored, with (6^^) oc 5^^, 



where z is an arbitrary direction. As described in Ref. [iTI, a variety of other condensates 
are, in general, also allowed. 

For large enough J, a second variety of superfluid is possible^, within which single 
bosons have not condensed. Instead, pairs of bosons condense, giving (fej^fej^) ^ 0, which 
does not break spin- rotation symmetry. For this state, the spin-singlet condensate (SSC), 




(2.2) 
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to be energetically favorable, J must be large enough to overcome the kinetic-energy cost 
associated with pairing. 

Referring back to the classes described in Section [T] and Fig. [T], we see that the transition 
from SSI into PC is an example of class A, in which the order parameter carries nonzero 
spin, while the transition from SSI into SSC is in class B, with a spinless order parameter. 
We will therefore mainly focus on these two phase transitions in the following. 



B. Strong-pairing limit 

To provide a concrete, if qualitative, guide to the phase structure of the particular model 
in Eq. (12.11) . we will implement a mean-field theory capable of describing the phases of 



interest. This will be similar to the approach of Ref. |20| for the spinless Bose- Hubbard 
model, where a mean-field is used to decouple the hopping term. 

Before describing this calculation, we will first use a very simple perturbative calculation 
to give an approximate criterion for condensation of boson pairs. In the limit of large J /U , 
an odd number of bosons on any site is strongly disfavored, and we can deal with a reduced 
Hilbert space of singlet pairs. 

The effective tunneling rate r for such pairs is given by r ~ t^/'^i, where vi = U + 2J 
is the energy of the intermediate state with a 'broken' pair. The effective repulsion of two 
pairs (ie, four bosons) on the same site is T ~ = 4t/. 

We therefore arrive at the simple criterion zt^ ^ U{U + 2 J) for the condensation of pairs, 
where z is coordination number of the lattice. This should be compared with the criterion 
zt > U + 2J for the condensation of single bosons^"^. These two simple results will be 
confirmed, and the numerical prefactors determined, by the mean-field analysis that follows 
(see Figure [1]). 



C. Quantum rotor operators 

To simplify the mean-field calculation somewhat, we will use S0(2) quantum rotor oper- 
ators hi and a^^ in place of boson operators in the mean-field calculation. These satisfy 

[oi^, fij] = Sijttif, (2.3) 

and 

K,ajj = 0. (2.4) 

This simplification, which automatically incorporates particle-hole symmetry, is convenient 
but inessential. The eigenvalues of hi are both positive and negative integers; they will be 
constrained to physical values by the potential. 
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The full Hamiltonian, from Eq. (12. ip . can be written a.s Ti. = V — T, where V is the 
on-site interaction and T is the kinetic energy term. In the rotor formalism, V is 

V = J2[U{n^-Ny + J\U\'] , (2.5) 

i 

where the term involving the chemical potential has been absorbed by making the spin- 
independent interaction explicitly symmetric about A^. For simplicity, we will take = 
2 in the following. In the rotor formalism, the angular momentum Lj is defined by its 
commutation relations with a^^, and the kinetic term is 

T = tY,{alaj, + h.c.) . (2.6) 

(id) 

First, consider the case when t = 0. Then the Hamiltonian is simply a sum of terms 
acting on a single site, containing only the commuting operators n and |Lp. The ground 
state on each site, which we label |2, L), is therefore an eigenstate both of n, with eigenvalue 
2, and of |Lp, with eigenvalue L{L + 1). For positive J, the ground state is a spin-singlet 
with L = and the lowest-lying 'charged' excitations^ are triplets with L = 1. 

On the other hand, for negative J, the maximal value of L is favored, which will lead 
(once a small intersite coupling is reinstated) to magnetic ordering. We therefore identify 
this latter case with the NI phase; the level crossing that occurs when J becomes negative 
will correspond in the thermodynamic limit to a first-order transition out of the orderless 
SSI. 

In the following we restrict to the case J > 0, to identify the phase boundaries into the 
SSC and PC phases. 



D. Mean-field Hamiltonian 



We will proceed by choosing a mean-field (variational) ansatz that incorporates the 
symmetry-breaking of the phases of interest. We choose to do so by defining a mean- field 
Hamiltonian TYmf, whose ground state will be taken as the variational ansatz. 

An appropriate mean-field Hamiltonian is 

^mf = — — — T$ , (2.7) 

where V is the same on-site interaction as in Eq. (12.51) and is the standard mean-field 
decoupling of the hopping term, generalized to the case with spin. 



^ Here, and in the following, we use the term 'charged' to refer to excitations that change the particle 
number, assigning to particles and holes charges of +1 and —1 respectively. Of course the bosons have no 
electric charge and feel no long-range forces. 
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where ip^ is a (c-number) constant vector, which will be used as a variational parameter. The 
remaining terms allow for the possibility of a spin-singlet condensate through the parameters 
^ and 

= Y.[^alal + ^*a,,ai,] , (2.9) 

i 

and 

T^ = Y. [^ala]^ + $*a,^a,^] , (2.10) 

where the sum is over nearest-neighbor pairs within the lattice. 

We now use the ground state of TYmf, which we denote |mf), as a variational ansatz and 
define 

E^f(V'^,^,<l>) = (mf|7^|mf) , (2.11) 

which should be minimized by varying the three parameters. If this minimum occurs for 
vanishing values of all three parameters, then |mf) breaks no symmetries and the SSI phase 
is favored. A nonzero value for ip^ at the minimum corresponds to PC, while vanishing ifj^, 
but nonzero values of \1/ and/or $ corresponds to SSC. 

Since Timf contains terms (within T$) that link adjacent sites, it cannot be straightfor- 
wardly diagonalized, as in the standard mean-field theory for the spinless Bose-Hubbard 
model. To find the phase boundaries, however, we need only terms up to quadratic order in 
the variational parameters, which can be found using perturbation theory. 



E. Variational wavefunction 



To order zero in if)^, \1/ and $, we require the ground state of V , Eq. (12. 5p . Assuming 
f/ > and J > 0, this is given by the simple product state 

|mf(°)) = ni2,0),. (2.12) 

i 

To first order, the ground state of Timf is 

Imfd)) = (-T^ + It* + |mfW) , (2.13) 

where Vi = U +2 J and f 2 = 4?7. (If rotors were not used in place of boson operators, a similar 
but somewhat more complicated expression would result.) All the physics incorporated in 
the mean-field ansatz is visible at this order: The first term allows for a condensate of single 
bosons, while the last two terms allow for a condensate of spin singlets. The third term is 
necessary to allow these singlet pairs to move around the lattice and make the SSC phase 
energetically favorable. 
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Computing to quadratic order in the variational parameters (which requires the 
expression for the perturbed states also to quadratic order) gives 
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(2.14) 



The transition to PC occurs when the top-left element in the matrix vanishes, while the 
transition to SSC occurs when the determinant of the remaining block vanishes. This gives 
the criteria 2zt > vi for PC and 27zt'^ > viV2 for SSC, in agreement with the simple 
considerations of Section IIIBi The phase boundaries were shown in Figure [H 

Note that it is in fact necessary to continue the expansion to fourth order to determine 
the direction of the vector ipfj_ when it is nonzero. This calculation has been performed in 
Ref. [2^, where the possibility of SSC was not incorporated and the PC phase was found to 
be favored, as expected. 

In principle it is also necessary to continue the expansion to higher order to investigate the 
competition between PC and SSC in the region where both are possible. Simple energetic 
considerations suggest, however, that condensation of single bosons in the PC phase will 
dominate over condensation of pairs, and this has been assumed in Figure [H 



III. SYMMETRIES AND PHASES 
A. Continuum action 

To describe the low-energy excitations of this model, we will derive a continuum field 
theory that captures the physics near to zero momentum. This assumes the absence of 
antiferromagnetic ordering of the spins, for example; it is chosen to be appropriate to the 
phases with which we are concerned. 

The action that results is in fact completely determined by the symmetries of the model, 
but a formal derivation is possible by analogy to the standard (spinless) Bose-Hubbard 
model^i2S. One first writes the partition function as a path integral and then decouples the 
hopping term using a site-dependent field ip^. Perturbation theory in tpf^ can then be used 
to eliminate all excitations above the ground state. 

The final form of the action has the same U(l) phase and SU(2) spin symmetry as the 
original Hamiltonian. Since the parameters have been chosen to give particle-hole symmetry, 
it can be written in a relativistic form: 

S^ = j d-^xdt {~4j,d'ij, + r4j^4j,)+S^^^ + ■■• . (3.1) 
Note that the action is completely relativistic and the derivative d acts in D = d + 1 
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dimensions: d"^ = — df. The quartic interaction contains two terms: 

= J d^xdt + . (3.2) 

The first term has 0(6) symmetry, while the second term, which vanishes if J = in 7i, 
breaks this down to SU(2). We will be interested in the case v < 0, for that yields a 
superfluid PC state. 

For higher spin F, the action has a similar form, with the field ip becoming a tensor of 
rank F. The quadratic part of the action is unchanged, but the quartic term now involves 
the F + 1 distinct scalar contractions of the field. 

B. Symmetries 

The action has full spatial symmetry, as well as the following (global) internal sym- 
metries. 

• Spin rotation (S*), under which ip^^ is a vector 

• Phase rotation ($): tp^ —>■ ■ip^e^'^ 

• Spin or phase inversion (/): ip^ — t/^^ 

• Time reversal (T): ip^ ^ ipfiy 9r ^ —dr 

These symmetries are not independent, and ground states breaking some of these symme- 
tries necessarily break others. Conversely, unbroken spin-rotation symmetry, for instance, 
implies unbroken spin-inversion symmetry, which we denote 

S ^ I. (3.3) 

Similarly, we have 

S ^ T (3.4) 

and 

$ ^ J. (3.5) 

Note that $ does not imply T, so it is possible to have broken time-reversal symmetry while 
maintaining phase-rotation symmetry. The identification of spin and phase inversion as the 
single operation / implies that it is impossible to break one without breaking the other—. 
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C. Observables 



To connect the predictions of the continuum field theory to physical observables, we must 
relate these to the field '^p^i■ 

Since angular momentum L is a pseudovector, symmetry requires 

Lp ~ kp^pijj^ip^ . (3.6) 

If time-reversal symmetry T is unbroken, we therefore have (L) = 0. (Note that unbroken 
/ is compatible with a nonzero (L), since L is a pseudovector.) 

In experiments with ultracold atoms, the population of each individual hyperfine state 
is also potentially measurable. We could consider, for instance, the operator h\j)iz, which 
counts the number of bosons in the /j, = z spin state (at site i). It is more convenient to 
define instead 

Qzz = blhz - lblb,p , (3.7) 

which measures the 'population imbalance' towards this state and has zero expectation value 
in a state without spin ordering. 

As suggested by the notation, Q^z is in fact a component of a (symmetric, traceless) 
second-rank tensor 

= \ [hlh. + hlhp) - ■ (3-8) 

In terms of the continuum fields, symmetry implies 

Qpu ~ ^ {'ippi^u + iJu^p) - ^s^uippipp ■ (3.9) 

D. Classification of phases 

We now list the phases described by the continuum theory of Section IIII At which allows 
for the breaking of spin and phase symmetries, but preserves the full spatial symmetry of 
the original lattice. The relevant connected correlation functions are the following: 



where (fo^f, ipo,i,iy cind ^I/o,i,f are scalar parameters. The unit vector rip is arbitrary and a^^^ 
and Up are defined as 

Opu = npUy - -5pu (3-11) 

and 

1 



^1. = ^{riip + m2p) , (3.12) 
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State 


Broken symmetries 


Nonzero parameters 


Spin-singlet insulator (SSI) 


None 


ipo 


Spin-singlet condensate (SSC) 


$ 




Nematic insulator (NI) 


S 




Strong-coupling pairing (SCP) 


s 


V^O, ^1, ^0, V'l 


Polar condensate (PC) 


/, $,5 


TpO, ipo, ^-1, ^'o, Ipl 


Ferromagnetic insulator (FI) 


S, T 


NI + 


Ferromagnetic SCP (FSCP) 


S, T 


SCP + Vf, 


Ferromagnetic condensate (FC) 


I, S, T 


PC + iff, il^i, ^{ 



TABLE I: The possible states described by the continuum theory of Section IIII A[ The notation 
for the symmetries is described in Section IIIIBI and the nomenclature for the states follows that 
used by Ref. |2l|. Note that the parameter ipQ breaks no symmetries and is therefore nonzero in 



every state. 



where ni and 112 are mutually orthogonal unit vectors satisfying rii x n2 = n. 

All of the states allowed by this symmetry analysis are summarized in Table [H Note that 
each state with S broken has a corresponding state with T also broken; these are states with 
at least one of (ff, ip{ and \l/f nonzero. Using Eq. (13. 6p . we see that they are ferromagnetically 
ordered and have (L) parallel to n. 

Table [T] is based on simple symmetry considerations and not every state listed will be 
possible in any particular physical realization. For example, the SCP phase appears, on the 
basis of energetic considerations, to be inevitably unfavorable compared to PC. Conversely, 
there are other possibilities for ordering that are not incorporated in Table [H such as breaking 
of spatial (lattice) symmetries ^. 

Since it is our purpose here to describe examples of the various classes identified in Section 
[U we will restrict our attention to a handful of phases and the transitions between them. 
Specifically, we will be interested in the three phases identified in Section [Tll SSI, SSC, and 
PC, allowing for the possibility of both spin and phase ordering. 



IV. PROPERTIES OF PHASES 

The various phases in Figure [U have different symmetries and low-energy excitations. We 
will now outline the properties of these phases, focusing in each case on the response to 
probes coupling to the angular momentum, L, and to Q^u, defined in Eq. (13. 8p . These are 



^ This is particularly important in Mott-insulating states of fcrniions, but may be less prominent in bosonic 
systems, where Bose enhancement encourages neighboring sites to have ferromagnetically aligned spins. 
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described by the correlation functions 

nL(x,t)~(7;L(x,t)-L(0,0)) (4.1) 

and 

nQ(x,t) ~ (7;g^,(x,t)g,,,(o,o)) , (4.2) 

where % denotes time-ordering. 

The exphcit calculations will be carried out two spatial dimensions, but most of the 
qualitative conclusions will also apply in three dimensions. 

A. Spin-singlet insulator 

The SSI phase is a 'featureless' insulator without spin or phase ordering. All quasiparticle 
excitations are gapped, ie, they occur at finite energy above the ground state. 

The phase can be further divided according to the lowest-energy 'charged' excitation. 
Throughout most of the phase, individual particle and hole excitations, described by the 
field V'/x, will have the smallest gap, but in a small region relatively close to the transition 
to SSC, bound singlet pairs will move to lower energy. (It is these excitations that condense 
across the transition to SSC, as described below, in Section IIVCI ) We will describe the 
former case here and return to the latter in Section IIVBI 

In the absence of any condensate, the appropriate action is simply that given in Eq. (13.11) . 
which we write, with r = A^, as 

u - , - , V - - 



Sssi= I d'^xdt 



'ip^ii-d'^ + A^)^^ + ^ ipi^i^f^ipiyipu + '^f^tpf^i'ui'u + ■■■ ■ (4.3) 



Particles and holes are described by the same field t/'^, with gap A and two distinct quartic 
interactions, with coefficients u and v. 

Perturbation theory in u and v, which we take to be on the same order, will be used to 
describe the low-energy properties of this phase. To do so, we first define the free propagator 
for the field V'/^: 

Because of relativistic invariance, the results are given below for the case k = 0; the corre- 
sponding expressions for nonzero momentum are given by the replacement u"^ — k"^ . 

1. Self energy 

First, we describe the self energy, for which the lowest-order 'tadpole' diagram is 



Sf = >. V >■ — (4.5) 
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FIG. 2: The spectral weight in SSI, calculated up to second order in the couplings u and v in 
d = 2. The delta-function peak at w = A (which has artificially been given a small but nonzero 
width) describes the stable particle excitation of the field 'i/'^i- For cj > 3A, there is continuum of 
three-particle excitations. 

(where the vertex represents some linear combination of u and v). This diagram does not 
depend on the frequency or momentum carried by the external line and so simply contributes 
a constant that renormalizes the gap A. 

The lowest-order diagram that depends on the external momentum is 



The physical interpretation of this diagram is as the decay of a particle, given sufficient 
energy, into a hole (described the top line, with the reversed propagation direction) and two 
particles. It is therefore clear that this diagram will make no contribution to the decay rate 
for a particle unless its energy exceeds 3A. 

This interpretation is clarified by using the spectral representation; the 'spectral weight' 
is given by 



where is the full propagator. This can be calculated numerically, and is shown in Figure 
[21 There is a delta-function peak'^ at ti; = A (a small width has been manually added to make 



We note here a subtlety regarding this plot and those in the remainder of Section Hvl The full series of 
diagrams, including those shown in Eq. (|4.5p and Eq. (|4.6p . cause a renormalization of the gap A (and 
hence movement of the features in p^) away from its bare value, ie, the value appearing expUcitly in the 
action. In this plot (and the ones that follow), A should be interpreted as meaning the renormalized value, 
rather than the bare value, and it is for this reason that the peak appears precisely sX uj ~ X. 




(4.6) 




(4.7) 
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it visible on the plot), and a three-body continuum, resulting from Eg , appears at lj = 3A. 
The spectral weight above this 3-body threshhold can be estimated from a non-relativistic 
theory of excitations above the gap, and yields p'^{uj) ~ (cj — 3XY~^ for u just above 3A. 



2. Spin response 

To describe the spin-response functions, we define 

^f.u,pa ~ {Tt t)M^, t)ijp{0, 0)^^,(0, 0)) , (4.8) 

in terms of which, 
and 

1 1 

Hq = - (n^^,^^ + ii^u,up) - -n^^,^^ . (4.10) 

The lowest-order diagram contributing to Ilfj,u,pa is given by the 'polarization bubble' 

no=----( )•■••, (4.11) 



which contains no interaction vertices. The dotted lines represent insertions of ip^ipi,; taking 
into account only this diagram gives spin-dependence Il^u,pcT oc 5^c,5yp. Performing the 
integration over the loop momentum in d = 2 gives the simple result 

no = ^log|^, (4.12) 

where uj is the frequency carried by the external lines. The imaginary part of this function 
is shown as the red line in Figure [3l and has a discontinuity at the threshold, as expected 
from Eq. 04.121) . In general the behavior just above threshold for the imaginary part is 
(cu-2A)('^-2)/2. 

Apart from insertions of the simple tadpole diagram into one of the propagators, 
which can be accounted for by a renormalization of the gap A, the only diagrams of first 
order in the couplings u and v have the form 





(4.13) 



Again, the vertex represents some linear combination of u and f , which we write u, and 
which depends on the particular response function of interest. Simple algebra in the spin 
indices gives u = | — f for IIl and u = | + f for IIq. In SSI, these two response functions 
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Im Hrpa 




u=l . 



1 . 5 



2.5 



3.5 



FIG. 3: The imaginary part of the spin-response function IIrpa in SSI in d = 2, calculated using 
the random-phase approximation. The red line, with zero coupling, gives the one-loop diagram 
Hq, while the other two lines give the results for nonzero coupling u. This coupling constant is a 
linear combination of u and v, the two couplings appearing in the action Sssi, Eq- (14. 3p . 



therefore have similar behavior, but are described by different combinations of the coupling 
constants. 

The diagram Hi is simply given by —u^q and in fact forms the second term of a geometric 
series. This set of diagrams can be summed, leading to the 'random-phase approximation' 
for the response function, 

Hrpa = . (4.14) 



1 + Mno 

In Figure [31 the imaginary part of this function is plotted in c? = 2, for three different values 
of the coupling u. It is nonzero only for u > 2A, the energy required to create a particle- 
hole pair. The discontinuity found earlier at threshold in Eq. fl4.12p is now suppressed 
logarithmically by the RPA corrections: and for small x = u — 2X > 0, the singularity is of 
the form 1/ logx. 



B. SSI near SSC 

As the spin-dependent interaction increases and the SSC phase is approached, a bound 
state composed of a singlet pair of bosons becomes energetically favorable. At the transition 
to SSC, the singlet pairs condense into a superfluid with no spin ordering. 

To describe the approach to this transition, we start with the action iSssi and introduce 
the field \1/ ~ 'ip^i'ip^i, by a Hubbard- Stratonovich decoupling of the quartic interaction v. The 
field \1/ describes the singlet pairs and will condense across the transition. It is described by 
the action 

S^= f d'^xdt [\d^\^ + r^\^\^ + ^\^\^+ . (4.15) 
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The full action has the form 

S'ssi =S^+S^ + ^j d'^xdt {^^f^^f, + ^^;.^;.) . (4.16) 

(The introduction of the field \1/ renormalizes the coupling constants within S^. Here and 
throughout, we will simplify the notation by retaining the same symbols for these renormal- 
ized quantities.) 



1. Self energy 



The diagrams shown in Section IIV A 11 coming from S^, will still contribute to the self 
energy near to SSC. As seen above, however, these diagrams are important only for u > 3A, 
whereas new diagrams coming from coupling to the pair field \1/ will contribute at lower 
frequencies. 

Using a double line for the propagator of the "if field, the first diagram is 




(4.17) 



in which a particle decays into a hole plus a pair. The vertices correspond to factors of 
g^. The threshold for this process is clearly A + A,^, where A is the gap to the particle and 
hole excitation as before, and A^ = y/r^ is the gap to the pair excitation. This excitation 
therefore becomes important at low energies when Aii, < 2A, which is simply the condition 
that a bound state exists below the two-particle continuum. 
The diagram can be evaluated in d = 2 to give 

^1 = ^ log T— r , (4-18 

and the corresponding spectral weight is shown in Figure HI the structure of this threshold 
singularity is similar to Eq. fl4.12p . and as was the case there, in general d we have a 
singularity ~ (cj — A — A^)*^'^"^)/^. A continuum of excitations appears for > A + A^, as 
expected. As the transition to SSC is approached, A^ becomes smaller and the edge of the 
continuum approaches the peak at u = A. The perturbation expansion used here breaks 
down as A,!, — *• and a more sophisticated RG calculation, described in Section IVIIt is 
required. 



2. Spin response 



Since \1/ is a spin singlet, it gives no direct contribution to the spin response, which is 
therefore given by Eq. (14. Sp . as before. The presence of the bound state, however, allows for 
new diagrams that contribute to the response function Ilf^u,pa- 
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Agi = . 05A 




0.8 1.2 1.4 1.6 1.8 2A 

FIG. 4: The spectral weight in SSI near the transition to SSC, calculated up to order g"^, plotted 
for three different values of A>ii, the gap to pair excitations. The coupling strength is g^, = 0.1. 
The peak at w = A (which is present for all values of A$ and has artificially been given a small but 
nonzero width) describes the stable particle and hole excitation of the field V'^i- For u > A + A$, 
there is continuum of excitations, corresponding physically to the conversion of a particle to a pair 
plus a hole. 

One such diagram is 




(4.19) 



which is the first term in the series 




(4.20) 



For X = io — 2X positive but very small, the pair propagator can be replaced by a constant. 
Summing over these diagrams will therefore lead to a similar 1 / log x singularity to that 
described above in Section [IV A 2[ This replacement is only valid for x < A^,. For u > 
2\ + X^^,, a 3-particle threshold singularity is found, similar to that discussed below Eq. (14.71) . 
Qualitatively different behavior will result as the transition to SSC is approached and A,j, — >• 
0. This will be addressed below in Section [VIII 



C. Spin-singlet condensate 



In the SSC phase, singlet pairs of bosons form a condensate, giving a superfluid with no 
spin ordering. This will only occur when the spin-dependent interaction is strong enough 
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to overcome the extra kinetic-energy cost of having bound singlet pairs. This requires 
sufficiently large J/U, as shown in Figured! 

The perturbation expansion used in Section IIVBI is not applicable here, and we must 
instead expand about the new ground state, with a condensed ^ field. (A similar approach 
can be used to describe the condensed phase of the spinless Bose gas2^.) 

We first write \1/ in terms of amplitude and phase as 

* = *oe'^' , (4.21) 

where \l/o and 9' are both real. For simplicity, we treat the amplitude of the \1/ field as a 
constant, ignoring the gapped amplitude modes. (This is appropriate sufficiently far from 
the transition to SSI, where the gap is large.) With this parametrization, iS^r, given by 
Eq. (I4.15p . can be rewritten as the action of a free, gapless field: 

S^ = j d'^xdt ^(9^)2 , (4.22) 

with the definition 9 = \/2^q9' . Physically, 9 is interpreted as the Goldstone mode corre- 
sponding to the broken phase symmetry in SSC 

In dealing with ip^, it is convenient to take out a factor of the condensate phase by writing 
V'/i ~ ^pL^^^' Then, since the condensate has broken phase-rotation invariance, we rewrite 
the field ip^ in terms of real and imaginary parts, 

(^, = i=(y.J + iv.J). (4.23) 
In terms of the new fields 9, ip^ and the action becomes'^ 

+ :y^yM - ^^H) ■ 90 + + ^i^i)id0r] . (4.24) 

Note that the two fields ip^ and (/:>^ remain gapped, but that their gaps are not the same. 
The lowest-energy 'charged' mode is (y?^, with gap A/ = \/r — 

1. Self energy 

As an example, we consider the lowest-order diagram that contributes to the decay rate 
for </)^, at energies well below that required to produce an excitation of the field ip^. This is 

^ It is straightforward to show that the Jacobian associated with the change of variables is equal to a 
constant. 
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-0.1 0.1 0.2 0.3 Ai 



FIG. 5: The spectral weight for the field in SSC, near the gap A/, calculated numerically 
using the diagram in Eq. (j4.25p . As in Figures [2] and HJ there is a delta-function peak at w = A/ 
(which has artificially been given a nonzero width), corresponding to the stable particle excitation. 
In this case, unlike in SSI, the continuum in the spectral weight occurs immediately above the 
peak. This is due to the (gapless) Goldstone mode resulting from the broken phase symmetry 
in SSC. The derivatives in the coupling between the Goldstone mode and the field in 5ssc 
strongly suppress the spectral weight as w — > A/ from above; in fact, ^ {uj — A/)^. 

given by 




(4.25) 



where the solid lines represent </?^ and the dashed lines 9. Each vertex represents a factor 
~ (pi ■ ^2)^0^; where pi and p2 are the momenta of the two 9 propagators, coming from 
the final term in iSssc- (AH the other interaction terms involve ^p^ and contribute to the 
decay rate only at higher energies. As in Section IIV A II there is also a lower-order tadpole 
diagram that does not contribute to the decay rate.) 

The diagram can be calculated numerically and the corresponding spectral weight is 
shown in Figure [5l As in SSI, there is a sharp peak (at uo = A/) corresponding to the stable 
gapped 'charged' mode, followed at higher energy by a continuum of excitations. In this 
case, however, the Goldstone mode 9 causes the continuum to begin precisely a.t u = Xj, 
albeit suppressed by a factor of {lj — Xj)^. 

This should be contrasted with the behavior at the transition itself, described below in 
Section rVlIi At the transition, the gapless modes are critical, rather than Goldstone modes, 
and their coupling is not suppressed by powers of the momentum. As a result, the spectral 
weight, calculated perturbatively, does not tend to zero as — A (see Section IVIIBp and 
a RG analysis shows that the sharp quasiparticle peak at = A is in fact replaced by a 
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weaker power- law singularity, reflecting the absence of quasiparticles at the critical point. 

2. Spin response 

The flrst contribution to the spin response resulting from the coupling of to 6 has 
three loops: 




(4.26) 



As for the self energy, the powers of momentum appearing in the interaction vertex strongly 
suppress the contribution at low energy, and this diagram will not modify the threshold 
singularity in the spin response. 



D. Polar condensate 

In the polar condensate phase, condenses, breaking both the spin- (S) and phase- 
rotation ($) symmetries. A different set of fields is therefore required to describe this phase, 
in which there are gapless excitations carrying both particle number and spin. A similar 
approach can be used to that described above for SSC, but since the calculations are rather 
involved, we will only give a brief outline. 

In terms of the general continuum theory of Section [III Al given in Eq. (13.11) . this phase 
corresponds to a condensate of ipfj,. This field is a complex vector and can take an arbitrary 
direction in space, which for simplicity we take as real and along the 2;-axis: (ipfj^) = <^q5^z- 
A convenient representation for the field is then to write 

V^^ = (/Poe^^e'°''^^e^^^-^ <5^, . (4.27) 

As in Section llVCt the amplitude modes of the condensed field, which are inessential to 
the physics, have been neglected. Summation over z,j G {x, y} is implied; Xi (^i ^"^^ Pj are 
real fields; and o",, and Tj are matrices acting in spin space: = ^ o o o^, Uy = ^ o o -i j , 

/ooi\ , /ooo\ 
r,. = and r,, = o o i . 

The physical interpretation of this parametrization is as follows: The factor e'-^ incorpo- 
rates the overall phase degree of freedom, similarly to e'^' in Eq. (I4.2ip . The matrices ai are 
the generators of real rotations in spin space and serve to rotate the axis along which the 
vector ijj^ is aligned. These two kinds of transformation correspond to symmetries of the 
action, and the energy is unchanged by uniform shifts in x, and ay. There are therefore 
three Goldstone modes in this phase, two of which carry spin. 

The remaining matrices, and Ty, generate complex transformations of the condensate 
wavefunction, and describe the angular momentum degree of freedom. The corresponding 



21 



modes, described by and Py, would be Goldstone modes in a system with full SU(3) 
symmetry and have gaps proportional to v, the coefficient appearing in Eq. fl3.2p . 
In this representation, the action can be rewritten 

5pc = I d^xdt !^^{dxf + lida^f + \ [mf - 4^/9'] - ■ {dx) + ■■■}, 

(4.28) 

where summation over i is again implied. [As in Eq. (14.221) . the fields have been rescaled by 
constant factors to give the coefficients of the kinetic terms their conventional values.] 

The angular momentum, given by Eq. (13. 6p . can be rewritten in terms of the fields defined 
in Eq. (I4.27p . giving for the transverse components 

Lx Py and Ly ^ (3^ , (4.29) 

and for the longitudinal component 

Lz ~ axPy - ay(3x . (4.30) 

The longitudinal and transverse components are given by different expressions because of 
the broken spin-rotation symmetry. These different expressions will lead to qualitatively 
different forms for the transverse and longitudinal spin responses. 

V. QUANTUM PHASE TRANSITIONS 

We now describe the behavior of the system across transitions out of the spin-singlet 
insulator (SSI) into the various phases described above. First we briefly address each of the 
transitions in turn, before describing them in more detail, in Sections IVIj IVIII and IVIIII 

First, consider the transition from SSI into the nematic insulator (NI). Both phases 
have {i})^) = 0, and off-diagonal elements of (^tl)^ipy) become nonzero across the transition. 
As suggested by the mean-fleld analysis Section III CI this is expected to be a first-order 
transition. This is confirmed by the presence of terms cubic in il)^il)y in the action S^, 
Eq. (13. ip . which are not forbidden by any symmetry (and are hidden in the ellipsis). We 
will not consider this transition further here. 

A. SSI to PC 

At the transition to PC, the field ip^ becomes critical; the appropriate field theory is 
therefore given by Eq. (14. 3p . Since the critical field ip^ carries spin, this transition is an 
example of class A identified above, in Section [H The field ip^j, is gapless and the spin 
response will be governed by the correlators of ipf_tipu, as in SSI. A RG analysis of this 
transition will be presented below in Section I VI I 
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B. SSI to SSC 



The critical field at the transition to SSC is the singlet pair introduced in Section [I VB[ 
Once has been isolated, ?/^^, which has no gapless excitations, can be safely integrated 
out. This leaves the field theory of a single complex scalar, with the same form as the action 
Sif, given in Eq. fl4.15p : 

= J d'^xdt (^|9^|2 + r^|^|2 + !^|^|4+ ...j . (5.1) 

This transition is therefore of the XY universality class, with upper critical dimension D = 
d + 1 = A. (The field \l/ has engineering dimension [\E'] = {d — l)/2, so the coupling has 
dimension [u^] = 3 — d.) 

Since the critical field \& is spinless, this phase transition falls within class B identified in 
Section [H While the action is sufficient to describe the critical properties of the ground 
state across the transition, we are primarily interested in excitations that carry spin, and the 
critical theory given by 5^ does not describe these. Instead, we must keep the singly-charged 
excitations given by ip^, and use the full action 5ggj, in Eq. fl4.16p . 

As in Section IIVBI the spin response is determined by the lowest-energy spin-carrying 
(but overall charge-neutral) excitations. These will normally be described by particle-hole 
pairs of ip^, but we also address the case in which there is a lower-energy bound state, which 
then governs the long-time response. Such a bound state is more likely to form in a region 
of the SSI-to-SSC phase boundary which is well away from the PC phase in Fig. [H 



1. Without bound state 



When there is no bound state, the spin response is determined by the correlators of the 
compound operator ip^ipu, as in Eq. (14.81) . The corresponding action is therefore given by 
5^si, in Eq. dUB]). 

Since ip^^ has only gapped excitations, while the field is now gapless, this can be 
simplified somewhat. The important excitations are those just above the gap A = ^/r, for 
which the dispersion can be replaced by a nonrelativistic form. We define particle and hole 
operators so that ^ + h^, giving an action Sm + iS(^ where 



2m„ 



A K 



2m^ 



(5.2) 



Using power counting (and taking [t] = [x] = — 1, since the critical theory Sxf, is relativis- 
tic), the engineering dimension of the kinetic-energy term is [1/m^] = — 1. The dispersion 
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is therefore irrelevant and, at least prima facie, the particles and holes can be treated as 
static impurities. We therefore take for the action 

S^,^ = j d'^x dt [p^ {idt + \)p^, + K {idt + \)h^ + {^PfJi^ + ^^m^m)] • (5-3) 

The scaling dimension of the coupling is [g^] = {3 — d)/2, so that it is relevant for 
d < 3. It is therefore relevant in two (spatial) dimensions and marginal in three, and we 
will consider both of these cases below. Any other interactions, including terms quartic in 
Pn and hfj^, are irrelevant. 

In Section IVIIt we will treat the case d = 2 using a renormalization-group (RG) analysis, 
and then, in Section IVIIBl return to the case d = 3, where straightforward perturbation 
theory is sufficient. 



2. With bound state 



If a bound state of and /ij, exists below the continuum, it will determine the lowest- 
energy spin response. At higher energies, the continuum of particle and hole excitations 
remains and will have the same effects as described above in Section IV B 1[ 

The particle and hole excitations and h^, carry spin 1, so the simplest spin-carrying 
bound states are a quintet with spin 2. (Bound states with antisymmetric relative spatial 
wavefunctions and net spin S* = 1 are also possible.) To describe these bound states, we 
introduce the field d^i,, a (symmetric, traceless) second-rank tensor: 

d^,u ~ -ipf^-ipu - ^Sf,uippipp . (5.4) 

By symmetry, angular momentum L couples to rf^i, by Lp ~ iepvpdp\d\y, and it is nec- 
essary to create a pair of d excitations to propagate the spin excitation. This will require 
more energy than an unbound particle-hole pair and so the response function IIl will still 
be given by the compound operator ippipu, as in Section IVB II 

On the other hand, it is clear from its definition, in Eq. (13.81) . that Q^^, couples directly 
to dp^. The response function IIq is therefore given by the two-point correlator of dp,^: 

Ugr^ {Ttdp,{^,t)dp,{0,0)) . (5.5) 

To construct the corresponding action, we notice first that there is no term such as ^"^dpp, 
because dpiy is traceless, so that the next order in the expansion must be taken. The action 
for dpi, is then given by 



Sq,d= d'^xdt 



V 



2 



dfiu idt - 7, ^ dpy + gddpudpy-^^! + 



The same power-counting argument as above shows that [l/m^] = —1 and the dispersion 
is irrelevant. The scaling dimension of the interaction gd is [gd\ = 2 — d, so that it becomes 
marginal ai d = 2 and is irrelevant for d = 3. This case therefore falls within class B2; we 
describe the properties of the critical region below, in Section IVIIII 
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VI. CRITICAL PROPERTIES: SSI/PC (CLASS A) 



The transition from SSI to PC is described by the action given in Eq. (13.11) . with two 
different quartic terms allowed by the U(l) ® 0(3) symmetry: 



— u — — V 

Sip = -i^f^d^i^^ + rip^ip^ + -%lj^ip^%ljyilj^ + -iJi,ip^,ilJyilJu . (6.1) 

To study the critical behavior of the model, a straightforward RG calculation can be 
performed, leading to the following beta functions to one loop: 

7 

13^ = -eu + -u^ + 2uv + 25^ (6.2) 
3 

= —ev + 3uv + -v^ (6.3) 

(where u and v are related to u and v by constant factors). These flow equations have no 
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stable fixed point for finite u and v. However, far more complete six-loop analyses in Ref. 
have shown that a stable fixed point does indeed exist with v < 0. At this fixed point, the 
full Green function for the field ipfj_ behaves like 

G^{p) r^p-^^"^^ , (6.4) 

where 77^ = to one-loop order, while the six-loop estimate is^i ^ 0.08. The polarization 
n^i/,pCT is determined by the response functions IIl and Uq defined in Section HVt and these 
have singularities of the form 

^L(g)~g-'+^^ Ugiq) ^ q~'+^<^ , (6.5) 

where the ?7l,q are exponents related to scaling dimensions of operators bilinear in ipfj^ at the 



fixed point of Ref. l2J. In particular, we have ?7l = f]H where the latter exponent is defined in 
Ref. l2J, and for which their estimate is ?7// ~ 2.70. For IIq we have rjQ = d + 3 — 2y^, where 



7/3 ^ 2.0 is the exponent listed in Table III of Ref. |25| for the collinear case with = 3. 



VII. CRITICAL PROPERTIES: SSI SSC (CLASS Bl) 

The critical theory for the phase transition between SSI and SSC is given by 

= j d'^xdr + r^\^\^ + ^1*1^) , (7.1) 

where we have rewritten the action in imaginary time. Since the field \i/ carries no spin, the 
lowest energy excitations are described by 

S^,^ = j d'^xdr [Pf, {dr + A) p^, + {dr + X)hf, + {^p^,h^ + ^^p^h^)] ■ (7.2) 
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For d < 3, the correlation functions of the particle and hole excitations can be found 
using a RG calculation. Since the present approach is slightly different from the standard 
RG, we perform the calculation using a cutoff in momentum space, which makes the logic 
involved more transparent, in Appendix |Xl Here we use dimensional regularization, which 
is the simplest approach from a calculational point of view. 

We define the (imaginary-time) free propagator for the \1/ field as 

G*(k,ia;) = - \ . (7.3) 

At the critical point, the renormalized mass of vanishes; in dimensional regularization, 
this occurs for = 0. For the and /i^ fields, the propagator is 

^^(i^) = — ^ , (7.4) 
—lu + A 

independent of the momentum. 

The renormalization of the terms in the action Si^, describing \& is identical to the standard 
analysis:— the presence of the gapped '0^ excitations cannot affect the critical behavior of 
the gapless \1/ field. The corresponding RG has a fixed point with uq, of order e = 3 — at 
which the scaling dimension of is [^] = 1 — e/2 + e^/100 + O (e^). 

To lowest order in the coupling (or, as will subsequently be shown to be equivalent, 
in an expansion in e), the only self-energy diagram for the particle field is as shown in 
Eq. dHT]): 

k, iu' 




= -"-^^--^--^^ (7.5) 



4/ r ^G^,{k,iu;')Gt{i{u;-u;')), (7.6) 




'k 



where 

d'^k 



oo 

d-1 



rid / dkk"-' (7.7) 



(for an isotropic integrand). 

There happen to be no diagrams giving a renormalization of the coupling at this order, 
but such diagrams appear in higher orders, as in Eq. 07.111) . 

Charge-neutral compound operators such as Lp and Q^^u, defined in Section UlI C[ can be 
written in the form T^^yp^hy, where T^y is a matrix of c- numbers. The critical exponents for 
these operators can then be found by considering the renormalization of the corresponding 
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insertion, given by 




j ^^G^{k,iu')Gt{iico-co'))Gt{i{co + co')) , (7.9) 

where all spin indices have been omitted in the latter expression. In accounting for these 
indices, it is important to note that the exchange of the pair interchanges the particle and 
hole lines and hence fi and u. This causes the results to be dependent on the symmetry of the 
matrix T^j,, leading to different scaling exponents for excitations with even (T^j, symmetric, 
eg, Qfiu) and odd spin (T^,^ antisymmetric, eg, Lp). 

To order e^, the diagrams that must be evaluated are, for the self energy: 




(7.10) 



and for the renormalization of g^: 




+ 




(7.11) 



(Note that the second diagram involves the four-point coupling of the \1/ field, indicated 
in the diagram by the empty circle.) The renormalization of the insertion is given by the 
following diagrams: 




+ 




+ 




(7.12) 



We will display the steps in the calculation of the results to order e and simply state the 
higher-order results afterwards. For the self energy, performing the integral over a;' using 
contour integration gives 



Sf(i^) 



9l 



1 

^2k 



1 



-ILU 



\ + k 



(7.13) 
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which leads to, defining z = —iu + A, 



sr(-) = . (7.14) 

plus terms that are finite as e ^ 0. 

To this order, the full propagator of the particle is then given by 

(G^)-! = 2 - (7.15) 

so that renormalizing the propagator (using minimal subtraction) at frequency z = fi gives 
for the wavefunction renormalization 

= ("7) 

Since there are no diagrams corresponding to renormalization of the coupling, we have 
Zg = 1, to this order. 

We now define the (dimensionless) renormalized coupling g^, given by 

= 27r~g^^^ . (7.18) 

In terms of this, we have Z^ = 1 — g'^/e. The beta function for the coupling is then given by 



dg^A 



so that the fixed point is at 



~9l = \ll- (7.20) 



Since the fixed point has g^ ~ e^/^, a perturbative expansion at this point is indeed equivalent 
to an expansion in e. 

The anomalous dimension of the particle (and hole) propagator is then given by 

V^ = f^:^\ogZ^, (7.21) 

0-9i) 

so that, to first order, 77^ = e/2 at the fixed point. 

The two-point Green function behaves, for u; > A, as 

G'^(cu) ~ (A-cu)"^+''* , (7.22) 

so that the corresponding spectral weight is given by 

p'^(cu) ~ (cu- A)-^+''* . (7.23) 
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FIG. 6: The spectral weight at the SSI-SSC transition in d = 2 spatial dimensions. The delta- 
function peak at w = A (present in the SSI and SSC phases) has been replaced by a continuum of 
excitations, with '-^ (uj — X)~^'^^^ . The numerical value rj^, = 0.91797 used in the plot results 
from a dimensional expansion in e = 3 — d, carried out to order e^, Eq. (j7.25p . and evaluated at 
e = 1. 

The relativistic invariance of the original theory allows these results to be extended to finite 
external momentum by the usual replacement u yZJ^^—W. 

The results at the next order in this expansion also involve diagrams renormalizing the 
coupling g^. The fixed point then occurs at 



i^V2 + (7 24) 

V15 looy ^ ' ' ^ ' 



and the anomalous dimension is 



Figure [6] shows the spectral weight p'^ for = 2 (e = 1), using the numerical value from 
Eq. (17.251) . The quasiparticle peak appearing on both sides of the critical points, Figures [2l 
m and 0, is replaced by an incoherent continuum of excitations. 

The scaling exponents for compound operators of the form T^^p^hy can be found in the 
same way, and the corresponding anomalous exponents are given by 

r7n±=/?;^log^, (7.26) 

for symmetric (+) and antisymmetric (— ) matrices T^y. The correlators of these opera- 
tors, n^, then yield the scaling dimensions of observables closely related to the 11 spin 
correlators defined in Section IIV A2[ This relationship will be discussed more explicitly in 
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Section lVlI A 2l below. where we will see that it is necessary to include the formally irrelevant 
momentum dispersion of the and particles to obtain the low momentum spin response 
functions of Section flV A 2[ From this argument we will obtain the general result that 

dim[n] = dim[n^] + d/2 , (7.27) 

where the corresponding 11 observables are taken on both sides of the equation. The corre- 
lation function is then given, for > A, by 

ni(2cu) ~ (A-cu)-^+''n± . (7.28) 

and the corresponding spectral density, 

A£(2cj) = lim Imn£(2cj + ie) , (7.29) 

is a delta-function at a; = A for r/n = 0, and otherwise behaves like 

A^{2ujf r^{uj- A)-^+^n± ^ (7 30) 

for uj just above the gap A. The physical spin correlation therefore has the spectral density, 
from Eq. (17:271) . 

A±(2cu) r^{uj- A)('^-2)/2+,n± _ (7 31) 

This result should be compared with the spectral density of Eq. fl4.12p . and the discussion 
below it, which holds both in the SSI and SSC phases. 

The perturbation calculation including diagrams up to two loops gives the exponent 

,n, = -.+ (^-^).^ + 0(.'). (7,32) 

for operators with even spin, such as Q^y, and r/n. = exactly for those with odd spin, 
such as the angular momentum Lp. The latter result is a consequence of the conservation 
of angular momentum. 

A. Effects of dispersion 

As argued in Section IV Bl the dispersion of the particle and hole are formally irrelevant 
in the RG calculation. Here we will present a few more details of the argument, and also 
show how the singularity in the physical momentum conserving spin response functions 11 
in Section flV A 21 are related to the local observables 11^ by Eq. fl7.27p . 
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1. Scaling form 

We first describe the standard scaling argument that suggests dispersion can be ignored. 
We concentrate here on the function A^y^p„^ which is sufficient to determine 11^,^ p^. via the 
Kramers-Kronig relations. (Similar considerations to the following apply to .) 

Consider the action of the RG transformation on the correlation function, exactly at the 
critical point. Let a = be the coefficient of k"^ in the quadratic part of the action of 
and hn and suppose that this is small (in some sense) but nonzero. Using real frequencies, 
with the shorthand 600 = u! — X, and suppressing the spin indices, scaling implies 

A{2{X + 5u), a) = fe^n^(2(A + b6uj), , (7.33) 

where Xn = 1 — ''7n and is the scaling dimension of the operator corresponding to disper- 
sion. By dimensional analysis, i/a = — 1 + 0(e), so the dispersion is irrelevant sufficiently 
close to dimension d = 3 including, we assume, at the physically important case of d = 2, 
e = 1. 

Using Eq. (17.33^ we can write the scaling form 

A(2{X + 6u), a) ~ 6uj^''^A{6u-y'^a) . (7.34) 

The power law in Eq. (17.301) corresponds to taking the limit a ^ in this expression, and — 
assuming analyticity at this point — is therefore valid when 6uj~'^°'a <^ 1. Since i/a < 0, we 
require that 6uj be sufficiently small. In other words, the dispersionless result for A[2ui) is 
appropriate for u just larger than A. 

2. Perturbation theory 

To provide a heuristic guess at the modifications caused by dispersion, we first consider 
the two-point correlator for the fields at finite momentum, ^^(k, cj). The simplest effect 
that can be expected is the replacement of the gap A by a momentum- dependent energy 
A + ak"^, corresponding to the free dispersion. 

This gives the slightly modified result 

G"^(k,cu) ~ (A + aA;'-u;)"'+''V (7.35) 

For fixed k, there remains a single peak, shifted from its k = position by a relatively small 
amount. 

To estimate the effect of dispersion on U^^^p^, consider the lowest-order diagram, evaluated 
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at external momentum equal to zero: 




n(°)(2M= ■••■( )■■■■ (7.36) 



/ Gt{K i(^ + u;'))Gt{-k, iiu - iv')) (7.37) 
^ ^ ^ (7.38) 



Note that Eq. (17.381) corresponds to replacing the dispersionless result, ~ (A — iuj)~^, by 
—iuj —iuj + a/c^, and summing over all momenta in the loop. It is clear that the same 
replacement can be made also in other higher order corrections with k replaced by other small 
momenta carried by the or /i^ quanta. We can heuristically account for such corrections 
in the above expression by writing (in real frequency) 

n(2cu) ^ j (A + «P _ ~^ + m gg^ 

"^dA;A;'^-i(A + aA;2-a;)"'+''" , (7.40) 



where a momentum cutoff A has been used. After taking the imaginary part, the integral 
can be performed to yield 

./ON f (cc^- A)('^-2)/2+r;n forA<cc;<A 

A{2uj) ~ < , ,1 , (7.41) 
^ ' \ (cj-A)-i+''n forcu>A ^ ' 

where A ~ aA^ is the bandwidth of the p^, excitations. Notice that the answer associated 
with the local spectral density, A^, appears at frequencies larger than the bandwidth, while 
the threshold singularity obeys Eq. (17.271) . 



B. At and above the upper critical dimension 

The critical results so far have been for d < 3, where they are controlled by non-zero 
fixed point value of g^. 

For (i = 3, the fixed point of the RG equations occurs for g'^ = 0, so perturbation theory 
in the coupling can be used to determine the structure of the Green function. Using the 
fully relativistic form of the action to calculate the lowest-order self-energy diagram gives, 
for > A^ + 

For comparison, the same calculation in two dimensions gives 

ImSt2 = — • . ^ (7.43) 
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FIG. 7: The spectral weight at the SSI-SSC transition in c? = 3 spatial dimensions. The 
delta-function peak at a; = A has been replaced by a continuum of excitations, with the mean-field 
exponent ~ (a; — A)~^. 

(ImE'^ = for u"^ < \^ + k'^ in both cases.) 

Note that Im S^^g tends to a constant as — + k"^ from above, in contrast to the 
cases considered in Section HVl above. The same is not true in three dimensions, and we have 

Im ~ - V)^T¥ , (7.44) 

and in general the threshold singularity is of the form (cu — \)^'^~'^\ For d < 3, this leading 
order estimate of the threshold singular is not correct, and it is necessary to resum higher 
order contribution by the RG, and was done in the previous subsection. For d > 3, higher 
order corrections are subdominant, and the leading order result here yields the correct 
singularity. Finally for c? = 3, we expect from Eq. fl7.19p that the coupling constant will 
acquire a logarithmic frequency dependence, g"^ ~ l/log(co' — A), and this will modify the 
above results by a logarithmic prefactor. The spectral weight for the particle and hole 
excitations in d = 3 (without the higher order logarithmic correction) is shown in Figure [7l 
As in Figure [6l the coherent quasiparticle peak is replaced by a continuum of excitations, but 
the exponent is given by its mean-field value: {u — A)"^. Note, however, that once the 

logarithmic correction has been included, the quasiparticle peak remains marginally stable 
m d = 3. The quasiparticle peak is well-defined for d> 3. 

Figure [8] shows Im Ti^ as a function of uj and fc, ioi d = 2 (upper plot) and d = 3 (lower 
plot). 
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VIII. CRITICAL PROPERTIES: SSI SSC (CLASS B2) 



The transition between SSI and SSC is, as in Section IVIIl described by the (imaginary- 
time) action 

S^ = j d'^xdr (^|9^|2 + r^|^|2 + !^|^|4j ^ (81) 

while the lowest-lying spin-carrying excitations now belong to the spin-2 field d^i,, with 
action ^ 

= y'd'^xdr d^u (^dr - ^ + ra^ df,u + gad^^df.u^'^ + ■ ■ ■ . (8.2) 

As described in Section IVB2t the dispersion is always irrelevant, while the coupling ga is 
marginal for d = 2 and irrelevant above. 

We therefore use perturbation theory in g^i to describe the response function Uq. Since 
Uq is related to the two-point correlator of dfj^^, by Eq. (15. 5p . we are interested in self-energy 
diagrams for d^jy. Diagrams such as 

= ^ ^ , (8.3) 

dfj,u d^y 

which do not depend on the external momentum or frequency, simply renormalize the value 
of Td and are of no interest. The lowest-order diagram that we consider is therefore 



d^u d^y 

Since \& is critical, rather than using the free propagator G* to describe this field, we use 
the full four-point (connected) correlation function, defined by 

n*(x, r) = (T. t(x, r)*(x, r)f (0, 0)^(0, 0)) - (0, 0)^(0, 0))' . (8.5) 

At the critical point, the Fourier transform of this correlation function has a power-law form, 
n*(k, iuj) = (7(vP+^)~^, where C is a constant. The exponent y is related to the scaling 
dimensions of the operator |\1/|2 at the critical point of iSij,; the latter scaling dimension is 
d + 1 — l/u (where u is the correlation length exponent of the SSI to SSC transition), and 
so y = 2/u — d—l. In d = 2, the current best valued is z/ ^ 0.6717, and so y ~ —0.0225. At 
leading order in e, a standard RG computation shows that y = |e + (9 (e^). (Note that using 
this one-loop result corresponds to summing the diagram shown for S^^"* as well as all ladder 
diagrams with increasing numbers of interactions between the two propagators.) 
The self-energy diagram S^^^ is then given by 

=9'J G',(i{uj + z.))n^(k, -iu) , (8.6) 
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where the free propagator for d^i, is given by 



Giiiu) = . (8.7) 

Using the exphcit form of 11* , we have 

r poo J 1 1 

The integrand is an analytic function apart from a pole at u = —uj — ir^ and a branch cut 
along the imaginary axis for u = ix, where \x\ > k. Deforming the contour along the branch 
cut in the upper half-plane gives 

4^^M = / fd^ . , . Im ,,, (8.9) 



TT 



k Jk 



X — iu! + (/C^ — X'^)yl'^ 



siCMny/2) f r^^^^.^^^^±^, ,3.10, 



Now consider the quasiparticle decay rate, found by analytically continuing iu; to real 
frequencies and taking the imaginary part, 

ImE<'>(r. + M = -,PM.y/2) ^^^^-^^0^^ . (8,11) 

where is the unit-step function. The frequency can now be isolated from the integral, 
leaving 

ImSW(r, + Su) = -gjC'/co'-y , (8.12) 

where Cy is a function of y only. 

Note that d > y, so that the decay rate, given by ImS^ evaluated at the quasiparticle 
energy, u = r^, vanishes. This is as expected, since the spin-2 field (i^j, is, by assumption, 
the lowest-lying spin-carrying excitation. We therefore conclude that, contrary to the case 
described above in Section IVIIl the quasiparticle peak survives, even in the presence of 
gapless critical excitations of \E'. 



IX. CONCLUSIONS 



This paper has described a variety of representative models of spin dynamics across 
the superfluid insulator transition of spinful bosons. The main classes can be discussed 
in the context of the simple mean field phase diagram in Fig. [1] of S" = 1 bosons, with 
the three phases SSI (spin-singlet insulator), SSC (spin-singlet condensate) and PC (polar 
condensate) . 

First, we presented the spin response spectral functions in the three phases. The SSI and 
SSC phases have a spin gap, and consequently, the associated spin spectral density A{uj) has 



36 



a non-zero value only above a threshold value given by the two ipfj_ particle continuum; the il^^ 
particle carries charge Q = 1 and spin F = 1. The nature of this threshold singularity was 
presented in Eq. (14 .121) and below, and below Eq. fl4.19p for the SSI phase. An alternative 
case, which could occur only well away from the PC phase, is that a ipf^ particle and an anti- 
particle form a Q = bound state with non-zero spin, d^j^; in this case, the spin spectral 
function in the SSI phase consists of a sharp quasiparticle delta function at the spin gap 
energy. Closely related results apply to the SSC phase, as discussed in Section ITV CI Finally, 
for the PC phase, we have a gapless spin excitation, leading to Goldstone spin responses 
noted in Section IIVDI 

Next, we described the quantum phase transitions between these phases under the fol- 
lowing classes: 

(A) This is the SSI to PC transition, associated with the condensation of the ipfj,. We found 
that the spin spectral density at the quantum critical point was determined by the scaling 
dimension of a composite spin operator which was bilinear in the ip^j^, as discussed in Sec- 
tion Ell 

(B) This was the transition from the SSI to the SSC phase, driven by the condensation of 
Q = 2, F = particle, This class had two subclasses. 

(Bl) The lowest non-zero spin excitation consisted of the two-particle continuum of the '0^ 
particle. This was the most novel case, and was discussed at length in Section IVIII Here we 
found interesting fluctuation corrections to the spin spectral density, characterized by the 
new 'impurity' exponent r]u in Eq. (17.321) . 

(B2) The spin response was associated with the sharp Q = 0, F = 1 quasiparticle d^^, which 
can be stable well away from the PC phase. The damping of this quasiparticle from the 
critical ^ fluctuations was found to be associated with a composite operator whose scaling 
dimension was obtained in Section IVIIII This damping led to powerlaw spectral absorption 
above the spin gap, but the quasiparticle peak survived even at the critical point. 

In all of the above phases and quantum transitions we also obtained results for the nature 
of the single-particle Green's function of the ipfj, particle. This is not directly associated 
with a spin oscillation, because the ipfj^ particle has a non-zero charge Q. However, this 
can be measured experimentally in 'photoemission' type experiments, such as microwave 
absorption, which involve ejection of one boson from the atom trap. 

We also did not consider the transition between the two superfiuids in Fig. [1], the SSC 
and PC. This transition is associated with spin rotation symmetry breaking, and so has a 
charge neutral, vector 0(3) order parameter. Consequently, the spin singularities can be 
mapped onto those of a relativistic 0(3) model, which were described in much detail in 
Ref. |28|. However, here we also have to worry about the coupling of the critical 0(3) modes 



to the gapless 'phonon' modes of the superfiuids. The nature of this coupling was discussed 
in a different context in Ref. |29|, and it was found that as long as the free energy exponent 
a < 0, the coupling between the Goldstone and critical modes was irrelevant. It is known 



that this is the case for the 0(3) model, and so the results of Ref. |28| apply unchanged to 
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the SSC to PC transition. 
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APPENDIX A: MOMENTUM CUTOFF RG 



In Section I VI It the scaling dimensions of the particle and hole excitations across the 
SSI-SSC transition were found using dimensional regularization. Here, we will peform the 
same calculation to one-loop order using a momentum cutoff. (The scaling dimension of the 
compound operator p^/i^ can be found by an analogous calculation.) 

Our approach will be to calculate the correlation functions of the gapped and /i^ 
excitations, evaluated for real frequencies just above the gap A. (Imaginary frequencies will 
be used as a formal device when calculating the diagrams, followed by analytic continuation.) 
We will find that there is a rescaling operation that is a symmetry of the theory and relates 
correlation functions evaluated at one frequency to those evaluated at another, as in a 
standard RG calculation. In this case, however, it is necessary to rescale relative to the gap 
energy A, rather than the zero of frequency. 



1. Self-energy renormalization 

As a result of particle conservation, there is no one-loop diagram contributing to the 
renormalization of the interaction vertex. 

The only one-loop diagram for the self energy of the particle (or hole) excitation is given 
in Eq. (US]): 

^tii^)=9i jj ^^GUKiu;')Gt{Ku;-u;')) , (Al) 

where 

(for an isotropic integrand), with A the cutoff. Since the dispersion of and is irrelevant, 
the diagram is calculated with the external momentum equal to zero. 
Working at criticality, where we set r,i, = 0, this gives 

after performing the integral over u' by contour integration. With the definition z = —iuj + X, 
we have 

EfM^^^/d*!^. (A4) 
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Using Dyson's equation, the inverse of the propagator is therefore 

Gtiiur = z-'-^ Tdk^ + O (4) . (A5) 



2 Jo k + z 

Taking the derivative with respect to A and expanding in powers of 2;/A gives 

A^ = _9|?fA-V»|^.A« + .... (A6) 

The first term is independent of u and so corresponds to a renormahzation of A, which is of 
no interest to us. The second term corresponds to wavefunction renormahzation and is the 
only relevant contribution from this diagram. 

Since there is no diagram giving a renormahzation of the coupling g^, a reduction in the 
cutoff from A to (1 — 6) A (with 6 infinitesimal) can be compensated by replacing the action 
5vi>,^ by 

5*,^ + SS^^^ = j d'^x dr I (1 + gl5) [p^ {d^ + X)p^, + (9, + A) h^] 

+ g^, {^p^,h^ + ^pf,h^) \ ■ (A7) 



To simplify this expression slightly, we have defined the dimensionless quantity^ g^ = 
^7^A-('^-3)/2^/lV2. 

2. Partition function 

This notion of 'compensating a reduction in the cutoff' can be made more precise by 
considering the partition function with discrete sources: 

ZA{Ji,'ki,Zi,g^) = / V^pV'^h exp-< / {p^zp^ + h^zh^) 

+ g^ j / [l'(ki - k2,i(t^i - u;2))p/x(ki,icui)/i^(k2,icu2) + c.c] 

+ '^^MfPA'lki, A - 2i) + h^{K A - Zi)] + c.c. i , (A8) 

from which correlation functions can be found by successive differentiation with respect to 
Ji and Jj. (We are concerned with ip^, so integration over with the appropriate measure, 
is implied.) The subscript A on the integral sign denotes that a cutoff A should be used. 



Note the similarity to the corresponding definition in Eq. (|7.18p . since fig ^ = 27r^. 
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Using this definition, Eq. (1A7|) can be written 



ZK{Ji,Zi,g^) = / P^pP^/i exp-<^ (1+^^(5) / {p^,zp^ + h^zh^,) 

J(1-5)A I J<^,k 

+ 94' jj {^Vi^hi. + c.c) + ^ Ji^[p^{\ - z.;) + h^{\- Zi)]+ c.c^ , (A9) 

which expresses the fact that the partition function, and hence all correlators, are unchanged 
by a shift in the cutoff and a compensating change in the action. A condensed notation has 
been used, where the momentum dependence is suppressed throughout. 

To bring this closer to the form of Eq. flASp . we rescale the fields and within the 
functional integral. By doing so, we can return the coefficient of the quadratic term to unity, 
giving 



ZK{Ji,Zi,g^) = I V'^pV'^h e^^-l / {p^,zp^ + h^zh^) 



'(1-5)A 

+ g^{l-gl6) II {^p^h^ + C.C.) 



+ X] (l - 1^5) Mp,{^ - + Ki^ - z^)] + c.c.j . (AlO) 

By comparison with Eq. (jASj), we can write 

ZAiJi^ki, Zi, g^) = Z(^i_s)a{{1 - gl5/2) Ji,^^, Zi, {1 - g^S) g^) . (All) 

3. Rescaling 

To return to the original theory, with cutoff A, we now perform a rescaling of all variables 
according to their engineering dimensions, with [r] = [x] = —1. Since we are working at the 
critical point of we have 

^(x/6, r/b) = 6^$(x, r) , (A12) 

where X is the scaling dimension of the field By dimensional analysis of Eq. (17. ip . the 
engineering dimension of \1/ is seen to be = {D — 2)/2 = {d — l)/2 and one would naively 
expect X = [^] = 1 — J. This expectation actually happens to be correct (to order e), since 
there is no wavefunction renormalization of S^s/ to one-loop order. 
Performing this rescaling leads to 

Za( Ji, K Zi, g^) = ZkAib^+^'/^Ju bki, bzi, b^-^'g^) , (A13) 

after making the substitutions \ — z) = b'^^'^^^ip{bk, A — bz) and X — icu' = {\ — iu)/b. 
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This can now be combined with Eq. (lAlip to give 



Za(J,, k,, Zi, g^) = Za((1 - ^>/2)(l + (1 + d/2)5).h, (1 + 5)k,, (1 + 5)zi, 

(1-^»(1 + (1-X)%^) . (A14) 

This gives a relationship between correlators in the same theory but at different frequencies. 
The fixed point of Eq. flA14p occurs when 

9i = [^-%?m + {^-m9i. (A15) 

so that 

gl = Vl^= (A16) 
which should be compared with Eq. fl7.20p . 

4. Renormalized propagator 

At the fixed point, Eq. (1A14I) becomes 

ZA(J„ki,2;„(7;) =Za((1+|/5)J„ (l + 5)k,, (1 + 5)2;,, , (A17) 

where y = 1 + d/2 — = 5/2 — 3e/4. Taking derivatives with respect to and 

gives 

('0/,(ki, A - Zi)tjj^{k2, A - Z2)) = 

(1 + 2y6){ij^{{l + 5)ki, A - (1 + S)z,)^,{{l + 5)k2, A - (1 + 5)^2)) • (A18) 

Using the conservation of frequency and momentum, we can define the propagator by 

(V^(ki, A - z,)M^2, X-Z2)) = {2nY+'6\k, - k2)6{z, - Z2)S^.G^{k,, X - z,) . (A19) 

Using this definition, Eq. (lAlSP becomes 



G^{k, A - 2) = (1 + y'6)G^{{l + 6)k, A - (1 + 6)z) , (A20) 

with y' = 2y-d-l = l- e/2. 

Restricting attention to k = gives 

G^{0, X-z) = {l + y'6)G^{0, A - (1 + 6)z) , (A21) 

which can be iterated to give 

G^{0,X- z) ^ z-^+'/^ . (A22) 
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Equivalently, after analytic continuation to real frequencies, we have 
which agrees with Eq. fl7.22p . 
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